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1 Introduction 



Yang and Lee pioneered a very interesting line of research in which one studies statistical 
mechanical models with the external magnetic field H generalized from real to complex 
values 0. Since the free energy of a spin system / = f(K,H) is a function of both the 
field H and the temperature T (or equivalently, K = j3J, where (3 = (/c B T) _1 and J is 
the spin-spin coupling), a related complexification is to generalize K from real to complex 
values. Both of these complexifications give one deeper insight into the properties of such 
models. Recently, there has been renewed interest in this subject for Ising models |2|]-[10| and 
Potts models ||11||-||13|| (earlier references can be found in these papers). As the comparison 
of complex-temperature phase diagrams for the 2D (isotropic, nearest-neighbor) Ising model 
with spin s > 1 versus s = 1/2 in Refs. [|j], |8||-||10| has shown, these phase diagrams differ 
considerably for different s even though all values of s are in the same universality class for 
the usual paramagnetic (PM) to ferromagnetic (FM) phase transition in this model. 

Another parameter of the Hamiltonian on which complex-temperature phase diagrams 
depend is the ratio of spin-spin couplings. Indeed, in cases such as the (spin 1/2, nearest- 
neighbor) Ising model on regular bipartite 2D lattices, where a variation in the ratio of 
spin-spin couplings along different lattice directions does not change the universality class 
of the PM to FM phase transition, this variation has a significant effect on the continuous 
locus of points where the free energy is non-analytic, which is denoted B: while B is a one- 
dimensional algebraic variety for the case of isotropic couplings, it becomes a two-dimensional 



variety for the non-isotropic case []14| . Recently, we have also shown, using exact results, 
that at complex-temperature singularities, the exponents describing the behavior of various 
thermodynamic functions depend, in general, on lattice type, which constitutes a violation 
of universality || . Moreover, we have found a number of violations of exponent relations at 
such singularities, such as 7 ^ 7' (for the susceptibility exponent) and a + 2(3 + 7 7^ 2 (as 
one approaches such a singularity from the PM phase) || . 

In the present paper, we explore further the extent of non-universal features of complex- 
temperature phase diagram of spin models, focussing on the effect of non-nearest-neighbor 
spin-spin couplings. Here, by "non-universal" , we mean features of the complex-temperature 
phase diagram that depend on a parameter in the Hamiltonian, where a variation of this 
parameter does not change the universality class of a phase transition at a given critical 
point. For our study, it will suffice to use a simple exactly solvable model, namely, the 
ID Ising model with nearest-neighbor (NN) and next-to-nearest-neighbor (NNN) spin-spin 
couplings, J nn and J nnn - This model (except for the special case Jnnn — — |<-^nn|/2) is critical 
at T = 0, so we specifically study the dependence of the complex-temperature phase diagram 
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on the ratio 

J nnn /, \ 

r = -j— (1) 

<Jnn 

for the range of r values where changes in r do not change the ground state of the model or the 



universality class of the critical point at T = 0. After some early papers [15, 16|, a detailed 



solution and discussion of this model was given by Stephenson in Ref. |17| and subsequent 
papers Some later papers include Refs. [[19], |20fl . All of these dealt with physical 

temperature; the model has not, to our knowledge, been studied for complex temperature. 
As our analysis of the complex-temperature phase diagram of this model will demonstrate, 
it illustrates very well how sensitive the CT phase diagram is to the presence of such NNN 
couplings. This example is also useful in giving one a qualitative idea, in a simple context, of 
what to expect concerning the effect of non-nearest-neighbor couplings in higher-dimensional 
spin models for which one does not have any exact solution. For these higher-dimensional 
models the addition of NNN couplings gives rise to quite complicated phase diagrams even 
for physical temperature |2l] . Moreover, although d = 1 is the lower critical dimensionality 



for the Ising model, and some features, such as the lack of a physical phase transition 
at finite temperature are qualitatively different from the behavior in d > dg_ — 1, past 
experience with 2 + e and 1 + e expansions [22] has shown that one can learn about spin 



models in intermediate dimensionalities by moving upward from di as well as downward from 
the upper critical dimensionality d u . Indeed, several intriguing connections were previously 
noted between the CT phase diagrams of the ID spin s Ising model and features of the 
same model on the square lattice [§, In passing, it may be noted that the ID NNN Ising 
model has also been used in a complementary manner in Ref. [p3fl , for a study of Yang-Lee 
(complex-field) zeros of the partition function for physical temperature. 



2 Ising Model and Notation 

In this section, we shall briefly give the relevant notation and review some basic properties 
of the model for physical temperature which will serve as a background for our new results 
on the complex-temperature properties. The ID (spin 1/2) Ising model with NN and NNN 
couplings is defined, for temperature T and external magnetic field if on a ID lattice, by 
the partition function Z = J2{a n } e _/3W with Hamiltonian 

T~t — —Jnn <J n cr n+l — Jnnn ° n® n+2 ~ H/,&n (2) 

n n n 

where a n = ±1 and (3 = (/c^T) -1 . Except where otherwise indicated, we take H = below. 
It is convenient to define 

K = (3J nn (3) 
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K' = (3J nnn (4) 

and h = f3H. Recall that on a bipartite lattice A, without loss of generality, one may take 
Jnn > (if J nn < initially, we can redefine o n — > —a n for n G A e , o n — > o n for n G A D , 
and J nn — > — J nn where A e and A Q denote the even and odd sublattices of A; the partition 
function Z is invariant under this mapping). Using this fact, we shall thus take J nn > 
henceforth. We define the ratio of couplings by eq. ([I]). We shall mainly consider the effect 
of NNN couplings with r > since (given that we take J nn > 0), these NNN couplings do not 
introduce any frustration or competition and do not change the universality class or ground 
state of the model. In contrast, given that we take J nn > 0, a NNN coupling with negative 
J nnn does introduce such competition and frustration and is not necessarily an irrelevant 
perturbation to the Hamiltonian. We shall also include some results for negative J nnn . It is 
convenient to define the Boltzmann weight variables z = e~ 2K ' , 

u = z 2 = e' iK ' (5) 

and z„ = e~ 2K , with 



K 

y 2 -4K 



u K = zl = e-- (6) 

For our study, it will be sufficient to consider the cases where (i) r = 1/p where p is an 
integer, or (ii) r is an integer, since these already amply demonstrate the sensitivity of the 
complex-temperature phase diagram to the value of r. In these cases, Z is a generalized 
polynomial (i.e. with positive and negative integer powers) in the respective Boltzmann 
weights (i) u and (ii) u K . Of course, one could also consider the case r = p/q where p ^ 1 
and g ^ 1 are relatively prime integers, but we shall not do this here. The (reduced, per 
site) free energy is defined as / = — (3F = lini7v- + oo A^ 1 In Z in the thermodynamic limit. 
We assume periodic boundary conditions and take the number of lattice sites A^ to be even 
in order to preserve the bipartite lattice structure on a finite lattice. 



3 Generalities and Complex- Temperature Phases 

For d — 1 dimension, it is straightforward to solve this model exactly, e.g., by transfer matrix 
methods. One has 

Z = Tr(T N ) =E A f ( 7 ) 

j 

where the Xj, j = 1, ...4 denote the eigenvalues of the transfer matrix T defined by T nn i = 
(£ n \exp(—/3E(£ n ,£ n '))\£ n i). It is natural to analyze the phase diagram in complex plane of 
the appropriate Boltzmann weight variables (such as u or u K for positive r). For physical 
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temperature, phase transitions are associated with degeneracy of leading eigenvalues. There 
is an obvious generalization of this to the case of complex temperature: in a given region of 
u or other Boltzmann weight variable, the eigenvalue of T which has maximal magnitude, 
X max gives the dominant contribution to Z and hence, in the thermodynamic limit, / re- 
ceives a contribution only from X max : f = ln(A max ). For complex K, f is, in general, also 
complex. The CT phase boundaries are determined by the degeneracy, in magnitude, of 
leading eigenvalues of T. As one moves from a region with one dominant eigenvalue X max to 
a region in which a different eigenvalue X' max dominates, there is a non-analyticity in / as it 
switches from / = ln(A max ) to / = \n(X' max ). The boundaries of these regions are defined by 
the degeneracy condition among dominant eigenvalues, |A ma:r | = |A^ a J. These form curves 
in the plane of the given Boltzmann weight variable. 

Of course, for physical temperature, a ID spin model with finite-range interactions has no 
non-analyticities for any (finite) value of K, so that, in particular, the ID NNN Ising model is 
analytic along the positive real axis in the complex u or u K plane and is only singular at T = 
0. In this context, we recall that the elements of the transfer matrix are non- negative (positive 
or zero) real functions of T for physical temperature, and the Perron-Frobenius theorem |24| 
guarantees that a (finite-dimensional, but not necessarily symmetric) square matrix with 
non-negative real entries has a real positive eigenvalue of greatest magnitude. This property 
underlies the absence of any non-analyticity and associated phase transition in a ID spin 
model with finite-range interactions. However, when one generalizes the temperature to 
complex values, the elements of the transfer matrix are not, in general, non- negative (they 
are complex), so that the premise of the Perron-Frobenius theorem is no longer satisfied, 
and, indeed, the maximal eigenvalue can switch as one varies K over complex values. 

Since the Xj are analytic functions of u, whence Xj(u*) = Xj(u)*, it follows that the 
solutions to the degeneracy equations defining the boundaries between different phases, 
| Aj| = \Xj\, are invariant under u u* . Hence the complex-temperature phase boundary 
B, or equivalently, the continuous locus of points where the free energy / is non-analytic, 
is invariant under u —* u* . The same applies for B in the u K plane, as discussed below. 
Although the model has a physical phase structure consisting only of the Z 2 -symmetric, dis- 
ordered phase, its complex-temperature phase diagram is nontrivial and exhibits a number 
of interesting features. 

Since the model has NNN couplings, one cannot define the transfer matrix as acting just 
between states consisting only of neighboring spins. The most compact way to define the 
transfer matrix is to use state vectors consisting of pairs of spins: v n = \a n , er n+ i). Then 

(v n \T\v n+ x) = (v n \e~ l3H \v n+1 ) = exp( — (a n a n+1 + a n+1 a n+2 ) + K'a n a n+2 + ha n+1 ) (8) 
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The factor of (1/2) is included because each interaction of spins within a given vector \v n ) 
is counted twice in the sum over n. Hence, with the basis vectors ordered as {\v n ) = 
I + +), I + -), I - +)> I - -)}, one has 0, |6 



T 



( e K+K '~ 


e -K>+ 



-K'+h 



-K+K'+i 







,-K+K'-h 


„-K'-h 





e -K'-h 



a K+K'-h 



\ 



(9) 



Note that T has zero matrix elements if the second spin in \v n ) has a value different from 
the the first spin in |f n +i), since these states overlap in this middle spin. Although T is not 
symmetric, the usual relation Z = Tr(T N ) = J2j , where the Aj's are the eigenvalues of 
T, still holds; this follows from (i) the theorem |24|] that an arbitrary complex i x i matrix 
can be put into upper triangular (u.t.) form by a unitary transformation V: VTV^ 1 = T u _ t ., 
such that diag{T u , t ) = {X x , ...X e }; and (ii) the identity Tr{T u N t ) = Tr(T N ) = £j =1 \f. It is 
convenient to define T = e -( K + K '+ h )j' so that T n = 1 and consider the eigenvalues of T. 
We shall take h = henceforth. 
T has the eigenvalues 



A 



i± 



A 



2± 



-K 



-K 



cosh K ± Vsinh 2 K + e~ 4K ' 



sinh K ± \/cosh 2 K - e~ 4K ' 



(10) 



(11) 



For physical temperature, Ai+ is the dominant eigenvalue, so the (reduced, per site) free 
energy is f — ln(Ai + ). In passing, we note that an equivalent method for solving the model 
in zero field is to re-express it formally in terms of a different theory with only NN couplings 
but a nonzero effective field [|T7|]. 

The internal or configurational energy (per site) U is 



U 



J nn sinh(lT) 



/sinh 2 (if) + e~ 4K ' 
Observe that 



'Jnn 



2e 



-4K' 



sinh 2 (i\T) + e~ iK ' + cosh(i\T) v /sinh 2 (i\:) + e 



-4K' 



(12) 



U(yJ n ni Jnnni P) U( Jnni Jnnni P) (13) 

which is an explicit illustration of the general fact noted above that we can, without loss of 
generality, take J nn > 0. The nature of the ground state (g.s.) depends on r |nj: if J nnn is 
positive or sufficiently weakly negative, the ground state is ferromagnetic: 



r > — 
2 



FM g.s. 



(14) 
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while for stronger negative J nnn it changes according to 

r<~\ (2,2) g.s. (15) 

where the (2, 2) g.s. refers to a spin configuration of the modulated form (-H H •••) 

Correspondingly, there is a non-analytic change in the ground state energy: 

U(T = 0) = -(J nn + J nnn ) for r>~ (16) 

whereas 

U(T = 0) = J nnn for r<- 1 - (17) 

Evidently (given that we take J nn > 0), negative values of J nnn give rise to competing 
interactions and frustration. Indeed, if J nnn is sufficiently negative that r < — 1/2, it changes 
the ground state of the model. In our study of the dependence of the complex-temperature 
phase diagram on the addition of irrelevant operators (i.e., irrelevant in the sense that they 
do not change the T = critical behavior), we therefore shall restrict to the case r > —1/2. 
However, since the range r < —1/2 is of interest in its own right, we shall also briefly 
digress to discuss this case further below. As Stephenson showed |17| 0, even in the range 



— 1/2 < r < 0, where the ground state still exhibits saturated FM long-range order, the 
NNN coupling has the interesting effect of giving rise to a "disorder temperature" Tp, where 
the correlation length has a local mininum; for T < To, the spin-spin correlation functions 
have a purely exponential asymptotic decay, while for T > To, their asymptotic decay is an 
exponential multiplied by an oscillatory factor. 

The T = criticality of the model is typical of a theory at its lower critical dimensionality, 
here d = 1. As T — > 0, the specific heat C has an essential zero given by 

C ~ 4Jfe B (l + 2r) 2 K 2 e- 2{1+2r)K as T -> for r > -- (18) 

and 

C ~ ^(i + 2r) 2 irV 1+2 ^ 

^ ^l( 1 -2\r\) 2 K 2 e~W r ^ K as T -> for r < -- (19) 
while for the borderline value r = — 1/2 one finds the proportionality 

C~k B AK 2 e~ 2K for r = -- (20) 
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For r > and for physical temperature, the spin-spin correlation function decays asymptot- 
ically like 

<^„>~(^) H (21) 



This is also true for —1/2 < r < OifT < Tu Hence, taking the T — > limit, 

one finds that the correlation length £, defined as usual by = — lim^oo n~ l In |(<Tocr n )|, 
diverges like 

e~ {l/2)e 2{l+2r)K as T^O (22) 
One also finds that for r > —1/2, the (zero-field) susceptibility \ diverges like 

X ~ (l/2)/T 1 e 2(1+2r)A ' as T -> (23) 

Consequently, for r > —1/2, 

C ~ , X ~ as T ^ , (24) 

independent of r in this range. Thus, for r > —1/2, the singularities in C and Xi expressed 
as functions of the correlation length £, are independent of r in this range, which shows that 
for r > —1/2, the NNN coupling is an irrelevant perturbation, and the model satisfies weak 
universality in the sense of Suzuki [p5], at the T = critical point. 

In order to investigate which eigenvalues are dominant in various complex-temperature 
phases, it is useful to express these as functions of the Boltzmann weight variables. For 
the case r — 1/p with integral p, since Z is a generalized polynomial in u, the CT phase 
diagram is well defined in the complex u plane. This follows since the CT zeros of Z may 
be unambiguously calculated in the u plane, and, in the thermodynamic limit, these merge 
to form the phase boundary B. In terms of the variable u, 

1 r 



A 1± = - 



1*2 ± J (1 — 1(2)2 + 4lt 1+ 2 



(25) 



1 -u* ± V(l + U2) 2 -4w 1+ i (26) 
Note that 

Ai± — > A 2 ± for y/u — > — s/u (27) 
If p is an odd integer, then ( p7|) implies that 



Ai±(w) = A2±(^*) = A2±(u)* for negative real u (28) 

respectively for the ± cases. (Here, we use the standard branch cut for y/u, along the 
negative real u axis.) 
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As background, we recall that for the case of nearest-neighbor couplings, B consists of 
the negative real axis in the u plane (e.g., Ref. |§). This is evident from the fact that for 
r = 0, the two nontrivial eigenvalues of T are Ai + = 1 + \fu and A2+ = 1 — \/u, which 
are equal in magnitude for negative real u. (The other two eigenvalues, Ai_ and A2- both 
vanish.) 

4 Case of r = 1/p for Positive Integer p 

We first consider the situation where the NNN coupling is of the same sign as, but weaker 
than, the NN coupling, i.e. < r < 1. For our purposes, it will suffice to deal with the case 
where r = 1/p with p a positive integer. There are two subcases: p even and p odd. For even 
p = 2£, the eigenvalues of the transfer matrix have the following Taylor series expansions 
about u = 0: 

A 1+ = l+M 1+ i + ... (29) 

\ 2+ = l-u 1+ 2 + ... (30) 

Ai_ = u§ + ... (31) 

A 2 _ = -u§ + ... (32) 

where ... denote higher-order terms in u. It follows that for even p, Ai+ is the dominant 
eigenvalue on the positive real u axis and hence also in the complex-temperature phase 
which includes this axis. Furthermore, on the negative real u axis in the vicinity of the 
origin, (i) if p = mod 4, i.e., £ is even, then A 2 + is the dominant eigenvalue, whereas (ii) if 
p = 2 mod 4, i.e., £ is odd, then Ai + is the dominant eigenvalue; in both cases, the respective 
eigenvalues are therefore also dominant in the CT phases which include this portion of the 
negative real u axis near the origin. For odd positive integral p, the A's have analogous series 
expansions in the z plane, Ai + = 1 + z 2+p + etc. Hence, in the u plane (with the usual 
+ sign taken for ^Ju if arg(u) = 0), Ai + is again dominant on the positive real u axis in the 
vicinity of the origin. 

Together with the theorem that a ID spin model with finite-range interactions has no 
non-analyticity for T > 0, i.e., along the positive u axis, it follows that for positive integral 
p, Ai + is the dominant eigenvalue on the entire positive u axis and hence the CT phase in 
the u plane which includes this axis and to which one can thus analytically continue from 
this axis. 
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We next prove a general theorem: For r = 1/p with p a positive integer, there are p + 2 
phase boundary curves emanating from the origin in the complex u plane, at the angles 

(2n+l)ir for 
2 + p 

Proof: To encompass the cases of both even and odd p, we use the Taylor series expansions 
in the z plane. From these, it follows that in the vicinity of z = 0, Ai + and A2+ alternate 
as the dominant eigenvalues. Now define polar coordinates according to z = p z e ldz , whence 
u = pe %e with p = p\ and 9 = 29 z ; then the degeneracy condition of leading eigenvalues, 
|Ai+| = IA2+I, reads |1 + z 2+p + ...\ = |1 — z 2+p + ...\ (where ... denote higher-order terms), the 
solution to which is cos((2 + p)B z ) = 0, i.e., 9 Z = (2n + 1)(tt/2)/(2 +p) for n = 0, ...,p + 1. 
This proves eq. (|33|). 

A related theorem is: For positive, odd p, the complex-temperature phase boundary B 
always contains the negative real u axis. To prove this, we again use the result that in the 
vicinity of z = 0, Ai+ and A2+ alternate as the dominant eigenvalues. We next observe that 
in the u plane, the degeneracy condition of leading eigenvalues, |Ai+| = IA2+I is automatically 
satisfied on the negative real u axis as a consequence of the symmetry condition (|27|) and 
relation ([28|). This completes the proof. 

Two further general theorems are the following: For positive r = 1/p with integer p, as 
one makes a half-circuit around the origin in the u plane, the dominant eigenvalues alternate 
between Ai+ and A2+. This is proved by noting first that from the Taylor series expansions 
above, these two eigenvalues are the dominant ones in the vicinity of the origin, and second, 
it is precisely their alternation as dominant eigenvalues which produces the phase boundaries 
emanating from the origin at the angles (|3~3"D and separating the different phases. Since the 
dominant eigenvalue at —9 is the same as that at 9, this theorem also completely determines 
the dominant eigenvalues on the rest of the full circle around the origin. 

By solving the degeneracy conditions of dominant eigenvalues, we have mapped out the 
complex-temperature phase diagrams. We consider odd values of p first and then even values. 
In Figs. [l]-|3| we show the results for r = 1, 1/3, and 1/5. 

For r = 1, we find that the complex-temperature phase diagram consists of three phases: 
(a) a region including the positive real u axis and extending outward to the circle at infinity, 
together with two complex-conjugate (c.c) phases (b), (b)*, located above and below the 
negative real u axis from the origin leftward to u = — 1. As follows from our general 
discussion above, Ai+ is the dominant eigenvalue in region (a), and, since A 14. and A 2 + 
alternate as dominant eigenvalues in the vicinity of u = 0, A2+ is dominant in regions (b) 
and (b)*. On curve separating region (a) from regions (b) and (b)*, |Ai+| = |A 2 +|. The CT 
phase boundary B includes multiple points at u = 0, where three curves meet, and a multiple 
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point at it = — 1, where four curves meet, with 2 different tangents, hence index 2. Here, we 
use the term "multiple point" and "index" in their technical algebraic geometry sense (see 
our previous discussions in Ref. ||). Thus, while in the model with only the NN spin-spin 
coupling, the CT phase boundary B is the negative real axis, the effect of adding a NNN 
spin-spin coupling with r = 1 is to produce two new complex-conjugate phases bounded by 
curves starting out from the origin and meeting at u = — 1. As must be true from general 
arguments, the theory is still analytic on the positive real u axis; the only changes are an 
increase in the number of CT phases elsewhere. 

The complex-temperature phase diagram for r = 1/3, Fig. |2], is progressively more 
complicated, consisting of five phases: (a) a region containing the positive real u axis; (b), 

(b) *, complex-conjugate phases whose borders are shaped somewhat like half-circles, adjacent 
to the negative u axis and including the interval — 1 < u < 0; and (c), (c)*, c.c. phases lying 
roughly concentrically outward from (b), (b)*, and including the interval of the negative real 
u axis —2.4 < u < —1. From our general discussion above, it follows that, in addition to 
region (a), Ai+ is dominant in regions (b) and (b)* while A2+ is dominant in regions (c) and 

(c) *. The CT phase boundary B contains multiple points at u = where five curves come 
together, at u = — 1 and u ~ —2.4, where in each case four curves come together with two 
different tangents, hence index 2. 

Finally, we show the complex-temperature phase diagram for r = 1/5 in Fig. |3|. The 
general features of this phase diagram follow from our previous discussion. A new aspect is 
that in addition to the part of B running along the negative real u axis, B also contains two 
complex-conjugate curves which extend to infinite distance from the origin in the "northeast" 
and "southeast" quadrants. 

We next show in Fig. f| a typical CT phase diagram in the u plane for an even value, 
p = 2, i.e., r = 1/2. In contrast to the diagrams with odd p, in those with even p, B does not 
contain the negative real axis. For p = 2, the CT phase diagram consists of four phases: (a) 
a region containing the positive real u axis and extending outwards to infinity; (b) a region 
including the interval — 1 < u < 0; and complex conjugate regions (c), (c)* above and below 
region (b). Our general discussion above determines the dominant eigenvalues in the various 
regions: Ai + in (a) and (b), and A2+ in regions (c) and (c)*. The CT phase boundary B 
involves multiple points at u = and u = — 1, each of index 2. The phase boundaries for 
even p may run to infinite distance from the origin. For example, we have also calculated the 
phase diagram for p = 4 case and find in this case that part of B consists of curves running 
to u = ±zoo. 
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5 Case of Positive Integer r 



We have also studied the situation where the NNN coupling is ferromagnetic and stronger 
than the NN coupling, i.e. r > 1; here, we focus on the case of positive integer r. For 
this case, the partition function is a generalized polynomial in u K . By re-expressing the 
eigenvalues Ai± and A 2 ± as functions of u K , i.e., 



A 



i± 



i+«l±V(i 



u 



l) 2 + 4ul 



+r 



\2± 



u 



I) 2 - Aul 



+r 



one sees that 



|Ai+(m, 



\^2+{ u k)\ , \\i-{u K )\ = \X 2 -(U I 



for real u K < 



(34) 
(35) 

(36) 



This implies that the CT phase boundary B contains the negative real u K . The expansions of 
these four eigenvalues around the origin of the u K plane follow directly from the expansions 
given in eqs. (|2T|)-(|3"2"D with the replacement u = u r . Hence, it is again true that as one 
traverses a half-circuit of the origin in the u K plane, the dominant eigenvalue along the 
positive real u K axis is A i+ and the dominant eigenvalues alternate between X 1+ and A 2+ . 
This also determines the dominant eigenvalues on the complex-conjugate half-circuit. These 
two results together imply that for arbitrary positive integral p, the CT phase boundary B 
always includes the negative real u K axis. 

Re-expressing the Taylor series expansions of the eigenvalues in terms of z K , the degener- 



acy condition for the leading eigenvalues, |Ai+| = 
where ... denote higher-order terms. Denoting u 

(2n+ 1)tt 



|A 2 +| reads |1 + z]^ 2r + . . . 



y l+2r 
-K 



+ ... 



Ad, 



the solution to this condition is 



l + 2r 



for 



n 



0,...2r 



(37) 



This proves that in the complex-temperature phase diagram in the u K plane, B contains 
1 + 2r curves emanating from the origin at the angles given in eq. (|37D . 

In Fig. ||we show the complex-temperature phase diagram for a typical case, r = 2. The 
dominant eigenvalues in the phases which are contiguous to the origin u K = are completely 
determined by our previous general results; starting from the phase containing the positive 
real u K axis and moving in the direction of increasing arg(ii K ), these alternate according to 
Ai+, A 2 +, and Ai + . For the remaining phase which is not contiguous to the origin, depending 
on one's choices of branch cuts connecting the branch points of the square roots in the 
eigenvalues, either A 2+ or A 2 _ is dominant. 
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Figure 5: Complex-temperature phase diagram in the u K plane for r = 2. 
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6 Negative r in the Interval — 1/2 < r < 



For the range — 1/2 < r < 0, notwithstanding the competition and frustration which the 
NNN interaction produces, it is still an irrelevant perturbation. As before, if 1/r = — p with 
positive integral p, then Z is a generalized polynomial in the variable u = e~ 4K ' and hence 
also in its inverse, w = 1/u. Since K' < 0, we shall use this inverse variable w for our 
analysis, in order to maintain the correspondence of zero temperature with the origin of the 
plots. The eigenvalues of the transfer matrix T in this variable can be obtained from eqs. 
(E3) and (EH): 

1 



Ai± = - 

\2± = \ 



1+w* ±J(l-W2) 2 + (38) 



W i±^{l + w2) 2 -A W 2- 1 (39) 



The borderline value p — 2, i.e., r = —1/2 at which the nature of the ground state changes, 
as discussed above, is evident in these eigenvalues since as p increases above p = 2, i.e., r 
decreases below r = —1/2, the eigenvalues cease to be finite at the origin w = because the 
last term in the square root becomes a negative power. By the same reasoning as before, if 
and only if p is an odd integer, B contains the negative real w axis. For integer p > 2, so 
that r > —1/2, which includes the region of interest here, these eigenvalues have the series 
expansions around w = 

A 1+ = 1 +W2- 1 + ... (40) 

A 2+ = l-«;|- 1 + ... (41) 

Ai_ = -w*- 1 + ... (42) 

A 1+ = wi- 1 + ... (43) 

where ... denote higher order terms. From this it follows that for our case of integer p > 2, 
Ai + is the dominant eigenvalue on the positive real w axis and in the CT phase which 
includes this axis. Other results are similar to those derived for positive r — 1/p above: on 
the negative real w axis in the vicinity of the origin, if p is even, then (i) if p = mod 4, 
then A2+ is the dominant eigenvalue, whereas (ii) if p = 2 mod 4, then Ai+ is dominant, and 
these respective eigenvalues are also dominant in the CT phase which includes this portion 
of the negative real w axis. 

In Fig. |6| and [7| we show our calculation of the complex-temperature phase diagram in 
the w plane for the values r = —1/4 and r = —1/3, respectively. For r = —1/4, besides (a) 
the wedge-shaped phase including the positive real w axis, where A i+ is dominant, there is a 
phase (b) which includes the interval — 1 < w < of the w axis, in which A2+ dominates, and 



17 




18 




19 



(c) a phase including the segment — oo < w < — 1 where A2- is dominant. In addition, there 
are two complex-conjugate phases (d) and (d)* where Ai_ is dominant; these have boundaries 
which cross each other at a multiple point of index 2 at w = — 1. For r = —1/3, besides the 
phase containing the positive real w axis, there are two pairs of complex-conjugate phases. 
As one moves from northeast to northwest, the eigenvalues which are dominant in these 
regions are A2+ and Ai_. 

For the present range —1/2 < r < 0, there is a finite physical disorder temperature To 



determined by the equation [17], [18 



cosh(if) = e~ 2K (44) 

The disorder temperature Tp decreases monotonically from Tjj = 00 at r = to Tjj = 
as r decreases to —1/2. In the context of the complex-temperature generalization of this 
model, we observe that, in addition to the physical solution of (pHf) for To, there are also 
complex-temperature solutions. In terms of the ratio r and the coupling K = Kr + iKi, the 
real and imaginary parts of eq. yield the respective equations 



cosh(^) cos(l£}) = e~ 2rKR cos(2rK 7 ) (45) 

sinh^) sin(fG) = -e~ 2rKR sm(2rK T ) (46) 
A more compact way of writing ( |44"|) is in terms of the Boltzmann variable z K : 

1 + z K = 2z r f> (47) 



Let us define K D = J nn / \k B T D ) and z K = e~ 2Ko . As an illustration, for r = —1/4, eq. (|47|) , 
expressed in terms of the variable u = w 1 ^ 2 , where w = e iK ' (whence, w(r = —1/4) = z^J 2 = 
e~ K ) is (uj — 1)(uj 3 +u) 2 +uj — I) = 0. (The trivial solution u = 1 corresponds to K = K' = in 
([44]) and is not of interest here.) The cubic factor has as roots the physical disorder solution 
uj d = 0.5437, i.e., wn = 0.2956 (K D = 1.219) and, in addition, the complex-temperature 
roots u= -0.7718 ±1.115z, i.e., w = -0.6478 ± 1.721* (K = -0.6094±1.931i). As noted in 
Ref . p| , there are an infinite number of complex K values corresponding to a given value of 
a Boltzmann weight variable, depending on one's choice of Riemann sheet in the evaluation 
of the logarithm; here, we list only one value of K for each w. These complex-temperature 
solutions of (0) lie in the phases (d) and (d)* in Fig. || From a similar analysis for r = —1/3, 
where w(r = —1/3) = e~ 4K ^ 3 , we find, besides the physical disorder point wo = 0.06694 
(K D = 2.028), also the two pairs of complex-temperature solutions w = —1.607 ± 1.539i 
(K = -0.5998 ± 1.783z) and w = 1.073 ± 1.366? (K = -0.4142 ± 0.6787i). One can see 
from Fig. |7| that these complex-temperature solutions of (0) lie in the interiors of four CT 
phases. 



20 



6.0 



4.0 



2.0 



3 0.0 
E 



-2.0 



-4.0 



-6.0 




-8.0 -6.0 -4.0 -2.0 0.0 2.0 4.0 6.0 8.0 

Re(u k ) 



Figure 8: Complex-temperature phase diagram in the u K plane for r = — 1. 

7 Negative r in the Interval — oo < r < — 1/2 

As discussed above, for r < —1/2, the NNN spin-spin coupling is so strong as to change 
the nature of the ground state from FM to the (2, 2) form. In Fig. [8] we show the complex- 
temperature phase diagram for a typical case, r = —1. Note that, in particular, by the 
same reasoning as for positive integer r (c.f. eq. (|36|) ), it follows that B always includes the 
negative real u K axis. 



8 Case r = -1/2 



For the borderline value r = —1/2, i.e., J nnn = — (1/2) J nn < 0, the competing preferences 
toward a ferromagnetic and (2,2) ground state are exactly balanced. Indeed, for r = —1/2, 
the model has nonzero ground state entropy, S(T = 0) = k B ln{(l/2)(l + V$)} 0, Hl> and 
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exponential asymptotic decay of (o"ocr n ) (modulated by an oscillatory factor) even at T = 
|T7|, [T8|1 . We find that if and only if r = —1/2, then the complex-temperature phase boundary 
B does not pass through the point T = 0, or equivalently, the origin in the w = 1/u = e 4K ' 
plane. This avoidance of the point w = by B shows the absence of criticality at T = 0. In 
Fig. P we present our calculation of the complex-temperature phase diagram for r = —1/2. 
One sees that B consists of two complex-conjugate curves which only intersect the real w 
axis at the point w = —1 (where they exhibit a multiple point of index 2). 



9 Potts Model 

Since the spin 1/2 Ising model is equivalent to the two-state Potts model, it is natural to 
extend the present study to include some remarks on how the complex-temperature phase 
diagram of the ID Potts model changes under the addition of a NNN coupling which is an 
irrelevant operator. Recall that in contrast to the 2D NN (spin 1/2) Ising model, no exact 
solution is known for general temperature of the 2D NN q state Potts model for q > 2 and 



hence the CT boundary B is not known even for this NN case (see, e.g., Refs. ]T2|, [13] and 
references therein). 

The zero-field ID g-state Potts model with NN and NNN interactions is defined by the 
partition function Zp = J2a n e _/Wp with 



Hp Jnn 3<j n <j n +i Jnnn 3<j„ a„+2 



where 8ij is the Kronecker delta and <r n G {l,...,q}. For our study, we shall consider ferro- 
magnetic couplings, J nn , J nnn > and define K, K', r as in eqs. ©-(0), and u p = e~ K . The 
ID NN model with J nn < involves finite ground state disorder, with ground state entropy 
So = ks ln(g — 1). This is similar to the situation on several higher- dimensional lattices; see 



Ref. |26| for the square lattice and Ref. for the honeycomb lattice. 

A notable feature of the CT phase diagram for the ID NN Potts model is its simplicity; 
in Ref. we found that for q > 3, it consists only of two phases, separated by the boundary 
B comprised of a circle 

- 1 + e iuJ 

u p = — , < lu < 2tt (49) 

Q ~ 2 

The solution of the ID NNN Potts model proceeds in the standard manner, via trans- 
fer matrix methods. As before, it is most efficient to use spin configuration vectors v n = 

|0n,0"n+l), SO that 

(v n \T v \v n+1 ) = (v n \e-P n \v n+1 ) = exp(y (S an ffn+1 + 5 an+1 ffn+2 ) + K'5 an an+2 J (50) 
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Thus the transfer matrix Tp is a q 2 x q 2 matrix. Here we consider the simplest case, q = 3. 
The resultant Tp is straightforwardly calculated from flSOl) . Defining Tp = e~( K+K ' +h ^Tp, we 
find for the characteristic polynomial of Tp: 

P{T P ;X) = (A + w p (1-w p ))(A 2 -(w p + w p + 1)A + w p (1 - w p )(1 + 2w p )) x 

x(A 3 + (u 2 p - l)X 2 + u 2 p (u p - l)(u p + 2)X + u 2 p (l-u p ) 2 (2u p + l)) 2 (51) 

The resultant eigenvalues of Tp are 

X = u p (u p -1) (52) 

1 r 



2 



1 + u p + u 2 ± Jl - 2u p - u 2 + 10« 3 p + < (53) 



together with three roots of the cubic factor in (0), each of which is a double root of 
P(Tp;X). We denote these as A3 a X^b, and A3 C and, since the expressions for these cubic 
roots are rather complicated, we omit listing them here. Aside from the polynomial Ao, the 
other eigenvalues have the following Taylor series expansions around u p = 0: 

X 1+ = l + 2u p + 0(u p ) (54) 

A x _ =u p +u 2 p + 0(u 3 p ) (55) 
X 3a = l-u 3 p + 0(u P ) (56) 

Xz b = -u p - l -u 2 p +0{u p ) (57) 
X3c = u p -~u 2 p + 0(u 3 p ) (58) 



In the absence of any NNN coupling, eq. ([49]) shows that B would be a circle of radius 
1 centered at u p = — 1. In Fig. |10| we show the complex-temperature phase diagram for 
r = 1. We find that the presence of the NNN interaction has a strong effect on this diagram. 
The CT phase boundary is much more complicated than just the unit circle centered at 
u p = —1. Rather than just two regions, as in the model with only NN spin-spin interactions, 
the complex-temperature phase diagram consists of nine phases. Three of these are (a) the 
region containing the positive real u p axis, where Ai + is dominant; (b) the region including 
the interval — 1 < u p < 0, in which A3 a is dominant; and (c) the region including the 
rest of the negative real axis, — oo < u p < —1, where Ao is dominant. The remaining six 
are comprised of three complex conjugate pairs. Starting from the northeast quadrant and 
moving to the northwest quadrant, the members of these pairs with Im(u p ) > are (d) a 
region with a wedge contiguous to the origin, where A3 Q is dominant; (e) a second region 
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Figure 10: Complex-temperature phase diagram of the ID q = 3 Potts model in the u p plane 
for r = 1. 

contiguous to the origin, where Ai + is dominant, and (f) at the same angle, ~ 27r/3, but 
farther out from the origin, a region where, depending on how the cuts linking the branch 
points of the cube roots are chosen, \ 3b or A 3c is dominant; the others are then the complex 
conjugates of these. The complex-temperature boundary B contains a multiple point at 
u p = 0, where six curves come together with three separate tangents (hence index 3), and 
two complex conjugate multiple points at (u p ) m = e 27U//3 and (u p )* m , where six curves meet 
in a tacnode, with three different tangents (see Ref. || for a discussion of tacnodes on CT 
boundaries £>). We note that the complex conjugate boundary curves separating phases (a) 
and (f) and (a) and (f)*, respectively, eventually head outward in northeast and southeast 
directions at larger distance \u p \ from the origin. 
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10 Conclusions 



In this work we have continued our exploration of the dependence of complex-temperature 
phase diagrams on details of the Hamiltonian, focussing on the effect of a next-nearest- 
neighbor spin-spin coupling in a simple exactly solvable model, the ID Ising model with 
nearest- and next-nearest-neighbor spin-spin couplings. Even for the range of values of 
r = Jnnn/Jnn where the NNN coupling is an irrelevant perturbation to the Hamiltonian at 
the T = critical point, we have shown that it has a considerable effect on the complex- 
temperature phase diagram. We have also presented some corresponding findings for the ID 
q = 3 Potts model. Our results further emphasize that, while complex-temperature phase 
diagrams and singularities give a deeper insight into the behavior of statistical mechanical 
models, they depend on details of the Hamiltonian, in contrast to the usual universality 
observed at physical critical points. 

This research was supported in part by the NSF grant PHY-93-09888. 
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